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We study the correlations in the one-dimensional extended quantum compass model in a trans- 
verse magnetic field. By exactly solving the Hamiltonian, we find that the quantum correlation of 
the ground state of one-dimensional quantum compass model is vanishing. We show that quan- 
tum discord can not only locate the quantum critical points, but also discern the orders of phase 
transitions. Furthermore, entanglement quantified by concurrence is also compared. 
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I. INTRODUCTION 

Recently, there has been a revived interest in the study 
of correlations in quantum systems. In the framework of 
quantum-classical dichotomy, the total correlations can 
be separated into a purely quantum part and a classical 
counterpart, and both quantum correlation and classi- 
cal correlation are able to be quantified respectively. A 
quantitative understanding of the different types of cor- 
relations might aid the application of quantum manipula- 
tion [T]. Therefore, distinguishing classical and quantum 
correlations in quantum systems is of both fundamental 
and practical importance. Paradoxically, entanglement 
was considered to be the most suitable manifestation 
of the quantum correlation and the main resource that 
speedup quantum computers over their classical counter- 
parts pi |3] . The studies of quantum correlations in ex- 
actly solvable models has a long tradition. For instance, 
the open-system dynamics of correlations in the presence 
of effect of environment is explored based on exact de- 
scription [3HS1- 

Let us consider a bipartite system described by the 
density operator pab shared by parts A and B, and if it 
can be written in the separable form pab = J2jPjPA^ 

p^g , such mixed bipartite state is termed disentangled. 
Though all the correlations therein are local, some sepa- 
rable quantum states still contain intrinsically quantum 
correlations 7J. In this sense, it seems entanglement is 
not always needed for quantum speed-ups [8 . Quan- 
tum discord (QD) was thus introduced to quantify non- 
classical correlations beyond entanglement paradigm in 
quantum states [5], and it has received an astonishingly 
amount of interest both theoretically and experimentally 

[Mill- 

In the field of quantum information, for a bipartite 
system pab, the total mutual information is the relative 
entropy between pab and pa® Pb, which corresponds 
to the minimal rate of randomness that is required to 
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completely erase all the correlations in pab, 

Hpab) ^ S{pa) + S{pb) - S{pab), (1) 

with von Neumann entropy S{p) = — Trplog(p). The 
quantum conditional entropy over a set of Von Neumann 
measurement is defined by 

S{pAB\Bk) -^^PhSipk), (2) 

k 

where the measurement-based conditional density oper- 
ator Pk associated with the measurement result k is 

Pk = —{Ia<E) Bk)pAB{lA ® Bk), (3) 
Pk 

in which I a is the identity operator of the subsystem A 
and pk= tr[{lA <8) Bk)pAB{lA ® Bk)]. Consequently, the 
associated quantum mutual information is given by 

I{pAB\Bk)^S{pA)-S{pAB\Bk). (4) 

The classical mutual correlation is defined as the superior 
of I{pAB\Bk) over all possible sets of one-dimensional 
positive-operator- valued measure (POVM) i?^, 

C(PAs) = sup /(pAsl^fc)- (5) 

The QD is then given by the difference of mutual infor- 
mation I(pab) and the classical correlation C{pab), 

D{pab)^I{pab)-C{pab)- (6) 

Since the minimization taken over POVMs is a notori- 
ous problem, so far only a few analytical results are ob- 
tained, including the Bell-diagonal states [El [16] , rank-2 
states 17 , and Gaussian states [IHl [TH]. Numerical ef- 
forts should be desired for general states. The QD can be 
shown to be asymmetric and nonnegative 120] , and it is 
invariant under local unitary transformations. The QD 
will vanish if and only if the state is classical, and it is 
implied that classical-only correlated quantum states are 
extremely rare 21j. 

The QD not only can discern classical and quantum 
correlations, but also can be used to establish relation 
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to quantum phase transitions (QPTs) of many-body sys- 
tems. Information on the locations and the orders of the 
QPTs can be obtained by consideration of the derivatives 
of the bipartite QD with respect to the coupUng param- 
eters [22], and has also been generalized to multipartite 
state |23l [24] . As we know, a QPT identifies any point 
of nonanalyticity in the ground-state (GS) energy of an 
infinite lattice system. The patterns in correlations of 
a many-body system suddenly change across the quan- 
tum critical point (QCP), and induce the non-analytic 
behavior of ground state |4'o)- The reduced density ma- 
trix (RDM) pij = Trij|*o)(^'o| is obtained by taking a 
partial trace over all degrees of freedom except particles 
i and j. The information of QPT is encoded in the non- 
analyticity of the matrix elements. With the QD being 
calculated from the RDM, one deduces that a discontinu- 
ity in the QD implies a first-order QPT, and a singularity 
in the derivative of the QD implies a second-order QPT 
[25| . Hence, the quantum correlation can serve as a hall- 
mark for the QPT in the many-body system. We would 
like to stress that the RDM comprises more accessible 
information than the information of QPT only. 

In this respect, we will take advantage of QD to study 
one-dimensional (ID) extended quantum compass model 
(EQCM) in the transverse magnetic field. The main rea- 
son for focusing on this model is that it not only allow 
us to study second-order phase transitions, but also al- 
low the investigation of first-order transitions. The rest 
of the paper is organized as follows. In Sec. [IT] we in- 
troduce the ID EQCM in the transverse magnetic field, 
and exploit the exact solutions. In Sec. Ill we calculate 
two-qubit QD of ID EQCM. We show that QD can not 
only diagnose various phase transitions, but also identify 
the character of the QPTs. Consequently, the zero tem- 
perature phase diagram of the model is mapped out. In 
Sec. |IV[ we recheck the entanglement in the ID EQCM 
in terms of concurrence. Sec. [V| finally contains the com- 
parison between the concurrence and the QD and a short 
summary. 



II. QUANTUM PHASE TRANSITION IN 
ONE-DIMENSIONAL EXTENDED QUANTUM 
COMPASS MODEL IN A TRANSVERSE FIELD 



z direction. N = 2N' is the number of the sites. We as- 
sume cyclic boundary conditions, i.e., the {N + l)th site 
is identified with the first site. The Hamiltonian describes 
a structure of two spins inside a unit cell. The orbital- 
orbital interactions depend strongly on the bond between 
two neighboring lattice sites. The Hamiltonian ([7| en- 
compasses two other well-known spin models: it turns 
into transverse Ising chain for Ji = Li, J2 = L2 = and 
the XY chain in a transverse field for Ji — Li,J2 — L2- 
The Hamiltonian ([7| can be exactly diagonalized by 
following the standard procedures. The Jordan- Wigner 
transformation maps explicitly between spin operators 
and spinless fermion operators by [5S1 [25] 



exp 



cr- = exp 



i-1 



i=l 



J-1 



CT- = 1 - 2c]c, 



'J - - (8) 

Next discrete Fourier transformation for plural spin sites 
is introduced by 



-1 = 77^^' 



k k 



= E^-"^^^' (9) 



with the discrete momentums as 

-(A^' - 1), -(A^' -3),...,N' -1. (10) 



Finally, the diagonalized form is achieved by a four- 
dimensional Bogoliubov transformation with two kind of 
quasiparticlcs 27., ,30,, 



1, 



(11) 



where optical spectra E^. — \/<^ + y^r and acoustic spec- 



tra El 



^. Here <j = \ak? + \Pk?. r 



(a*./3fc +afc/3*)^ +4|afc|^/i^ a,. = ( Ji + J2) + (Li +L2)e'*, 
and /3fc = (Ji — J2) — (Li — £2)6'*^. The ground state Eq 
is obtained, 



(12) 



The Hamiltonian of ID EQCM in an external trans- 
verse magnetic field is given by [26l [21] 

N' 



(7) 



where erf (a = x, y, z) is the Pauli operator at site z, 
Ji and J2 ( Li and L2) are the strength of the nearest- 
neighbor interaction on the odd (even) bond, and h char- 
acterizes the intensity of the external field applied in the 



It is easy to find that t he energy gap of ac oustic branch 
will close when h = 2^/ ( Ji ± L2)(>^2 i Li), and the non- 
analyticities of the ground state determine QCPs. In the 
absence of the transverse magnetic field, i.e., h = 0, the 
critical lines correspond to Ji = ±£2 and J2 — ±ii, 
respectively, which confirms the conclusion in Ref. [31) . 



III. TWO-QUBIT QUANTUM DISCORD OF 
EXTENDED QUANTUM COMPASS MODEL 

The QD is explored from the two-qubit RDM. In the 
representation spanned by the two-qubit product states 
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{|0)a®|0)b, \0)a<»\1)b, \1)a'S)\Q)b, where 
|0) (|1)) denotes spin up (down) state, the two-site den- 
sity matrix can be expressed as. 



3 

El „a \ „a a' 



(13) 



where erf are Pauh matrices erf, a\ and af for a = 1 to 
3, and 2 by 2 unit matrix for a=0. The Hamiltonian has 
Zi symmetry, namely, the invariance under parity trans- 
formation P = (S)iO'l , and then correlation functions such 
as {(jfTj) {a — x,y and b = 0,z) simultaneously vanish. 
Also, {(jfcTj) ((o-f crj)) is zero due to the imaginary char- 
acter of crj (cr^). Therefore, the two-qubit density matrix 
reduces to a X-state , 




(14) 



with 



7.± = -(l±2(a,f) + (afa|)), 
z± = i(«aj)±(a,rap), 

Wi = J(l- (crfcr|)). 



(15) 
(16) 

(17) 



The density matrix of single qubit is easily obtained by 
a partial trace over one of the two qubits. 



(18) 



Thus the total correlation is quantified by the quantum 
mutual information as 



HPij) = S{p^) + S{p-i) - S{p^J), 



(19) 



with S{p.)=Sip,)= -E™=o {[1 + 

(-l)™«)]/2}log2{[l + (-l)"(a,f)]/2} and 

S{Pij)^- EL=0 ^™ l0g2 Cm - Eri=0 C" ^°S2 ^n, whcrC ^„ 

= [l + (afa|) + (-l)'"^((afaj) - (araj))2 +4«)2]/4 

and = «^|) + (-l)"((afap + (afap)]/4. 

Since part B contains a single qubit, we can compute 
the classical correlation by extremizing Eq. ([5]) over a 
complete set of orthogonal projectors {Bk = |0k)(6k|, 
K =11, -L}, where Q\\ = cos(6l/2)|0>i3 + e*'^ sin(e'/2)|l)B 
and e_L = e-*'^sin(e'/2)|0)i3 - cos(e'/2)|l)B with 
{) < 9 < TT and < < 27r. Extensive numerical analysis 
implies that the extremization is mostly achieved at 
9 = TT /2, (f) ~ Q [32 [33], which is also confirmed in 
our numerical optimization. Accordingly, the classical 
correlation is expressed as 



where ffbin(p) = -plog(p) - (1 - p) log(l - p) is 
the binary entropy, pi = (1-1- (crf))/2 and p2 = 

,z\2\ /2. Thus the 



1-f ./fmax(|(afaj)|,|(afaj) 



quantum correlation is simply given by 
D{Pij) = I{p^j) ~ C{p,j). 



(21) 



The Eq.pO]) has been verified in Werner state [H], ID 
Ising model, XY model [34], and XXZ model [22l |32], 
even for finite temperature |34H36j . 

First, we consider the QPT in ID XY model, i.e., 
Ji = Li, J2 ^ L2. 7 = (Ji - J2)/{Ji + J2) denotes 
the anisotropy of the coupling. The phase diagram is 
shown in Fig. [T] which describes two distinct QPTs. The 
transverse magnetic field h drives a second-order transi- 
tion from an antiferromagnetic ordered phase to a para- 
magnetic quantum-disordered phase, and the anisotropy 
critical line 7 = is the boundary between a Neel state 
along the x direction and a Neel state along the y di- 
rection [37l [38]. The first derivative of the QD for the 
nearest-neighboring odd bond is displayed in Fig. [l] It 
is evident there are rapid changes of QD in the critical re- 
gions. The changing rate of the QD is most pronounced 
for critical line 7 = 0. The diverging first derivative 
of QD suggests both transitions belong to second-order 
phase transitions. Finite size scaling shows that the peak 
will not increase with respect to lattice size. 

Next, we consider the QPT in ID compass model, i.e., 
Ji = L2 = 0. The ground state of finite-size system is 
2^ fold degenerate [59]. In the absence of magnetic 
field ft, = 0, we note that there is only classical corre- 
lation cre ated between nearest neighbors. In this case. 



2J2LiCosfc, El = 0. The only cor- 
relation functions surviving are C|j_j^ 2i = {'^2i-i'^2i) 
^21 2i+i = (^2i^2i+i)- As a consequence, the density ma- 
trix becomes densities diagonal in the orthogonal prod- 
uct bases. In other words, they are equivalent up to local 
unitary operations to Bell-diagonal states [HTj . 



P2i-l,2i 



/I 








^2i 





1 


^2i-l,2i 








nV 

'-'21-1, 2i 
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-1,2. 





1 



2i \ 
I 

(22) 



C(pij) = -ffbm(Pl) - Hun{P2), 



(20) 



where are four Bell states with |\E'i) = (|00) + 

|11))/V2, 1*2) = (101) + |10))/y2, IM/3) = (101) - 
|10))/y2, 1*4) = (|00) - \n))/V2, and Ai - A3 = 
(1 - C|.-i.2.)/8, A2 - A4 - (1 + C|,_i,2.)/8. We can 
find that QD vanishes in case of maximum Ai is less 
than 0.5 [HHlj, that is, -D(p2i-i,2i)=0. Similarly, the 
quantum correlation between two qubits on even bonds 
is also equal to zero, i.e., -D(p2i,2i-i-i)=0. In a sense, 
ID compass model behaves as classical system. The 
null QD is induced by the macroscopic degeneracy of 
the ground state dut to the peculiar symmetry of the 



4 




Critical Ising 



-1 



r 




FIG. 1: (Color online) Top: phase diagram for the anisotropic 
XY spin-1/2 chain in a transverse field at zero temperature. 
The heavy lines represent second order phase transitions. The 
horizontal line will be referred to as the anisotropic transition 
and the vertical line as the Ising transition. PM denotes a 
paramagnetic phase and AFM^, (AFMy) denotes antiferro- 
magnetic phase along the x (y) direction. The XX model 
obtained by setting 7 = displays a critical line for magnetic 
field h £ [0,2]. The Ising model obtained for 7=1 exhibits 
a critical point at /i = 2. Bottom: the first derivative of the 
QD between nearest neighbors in the 7 — ft plane. 



Hamiltonian [12] ■ The similar strategy applies for the 
case J2 = 0, L2 = 0, h = 0. As a consequence, 
the energy spectra are greatly simplified as = 2Ji 
and = 2Li. Then, the correlation functions are 
found to be (crf,_i0-|i) = (cr2zcrfi+i) = and (crf,_if7f,) = 
(crfjtjfj^i) = —1. In such case, D{pi^i+i) — 0. We have 
verified that if a bipartite quantum state is a product 
state, i.e., pab = Pa ® Pb, the state has no quantum 
correlations, not vice versa. 

It has been shown that zero QD between a quantum 
system and its environment is necessary and sufficient for 
describing the evolution of the system through a com- 
pletely positive map [43l |44]. In addition, a quantum 
state can be locally broadcast, i.e., of locally sharing 
preestablished correlations, if and only if it has zero QD 
[13 m] . These purely classical states are rather unsteady, 
and a generic arbitrarily small perturbation will make 
them become nonclassical [H]. 

In the following, we will focus on the case of J2 > 0, 




FIG. 2: (Color online) Top: QD of nearest-neighbor sites on 
odd bond as functions of Ji and J2 with parameters Li = 
1,1/2 = 0, ft = 0,iV = 1024. Bottom: the derivative of the 
QD in the Ji — J2 plane. 



Li — I. Minus sign of J2 is irrelevant because of unitary 
equivalence. In absence of external magnetic field, i.e., 
h = 0, there is a first-order QPT at Ji = which sep- 
arates the phases with (cfj-icrfj = 1 from those with 
{o'2i-i<^2i) — There is also a second-order QPT at 
J2 — 1 separating the phases with (trf ) — from those 
with (af) 7^ 47 . Two critical lines where the energy 
gaps vanish, separate four gapped phases in the parame- 
ter space, as is shown in Fig. [2j Using Eqs. (191 - (21 1, 
the correlations can be evaluated in detail. We observe 
that QD indeed vanishes for Ji = 0. The carmine con- 



tour line (very large values beyond the scope of colorbar) 
implies a first-order critical line while yellow contour line 
implies a second-order critical line. The occurrence of 
discontinuous dips at Ji = implies a first-order critical 
line and a singularity of the first-order derivative of the 
QD reveals a second-order phase transition. To interpret 
the relation between the QD and QPTs more clearly, let 
us examine QPTs along three paths which start at the 
point Ji=l, J2=0, where Hamiltonian reduces to that of 
the ID quantmn Ising model. See Fig. [3] for an illustra- 
tion. 

In Fig. |4] we plot the QD verse Ji along the three 
paths shown in Fig. [3j It is clear that QD exhibits a 
sudden drop at Ji — 0, suggesting a first-order QPT. For 
paths J2 = (1 — Ji)/2 and J2 ~ 2(1 — Ji), a diverging first 
derivative of QD signals a second-order QPT. As for the 
multicritical point on the line J2 = 1 — Ji, a first-order 
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FIG. 3: (Color online) Phase diagram of the ID quantum 
compass model in the absence of magnetic field. A first-order 
QPT occurs at critical line Ji — and a second-order QPT 
at J2 = 1 113. We plot QD in Fig. Q along three paths (a), 
(b), and (c) indicated by the dashed lines, which start at the 
point Ji — l, J2~0, where the Hamiltonian ([7| reduces to that 
of the one-dimensional quantum Ising model without external 
magnetic field. 



phase transition dominates indicated by a discontinuity 
of the QD. Note that the classial correlation C{pij) = 1 
at J I ^0, and also has an abrupt decrease at Ji =0. 

As mentioned above, the zero discord state is incred- 
ibly fragile, and finite perturbation can create QD. For 
example, the magnetic field will lift the high degeneracy 
of the ID compass model. Therefore, we are especially 
interested in the role of the QD in detecting QPTs in 
ID extended compass model under the magnetic field. 
In Fig. |5] we plot the QD of nearest-neighbor qubits as 
the function of Ji and J2 with external magnetic field. 
We present the contour map of the first derivative of the 
QD. With the increase of magnetic field h, we find that 
the critical lines deviate from those for h — and the 
second-order critical line develops into hyperbolas. 

Furthermore, we plot QD as a function of Ji along the 
path J2 = Ji and its corresponding first-derivative for 
/i = 2 in Fig. [6) The divergent peaks s how clea rly that 
there will be QPTs at Jic = (Li ± ^jLl + h'^)/2 and 
(— Li ± Y/Lf^^7i^)/2, and these QPTs belong to second 
order. Hence, the phase diagram of ID EQCM in the 
transverse magnetic field is sketched in Fig. [Tj which is 
identical with that analyzed by correlation functions |30| . 
It is worthy noting that there are only continuous phase 
transitions once applying the magnetic field. 
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FIG. 4: (Color online) (a) The total mutual correlation, the 
classical correlation and quantum discord along the paths (a) 
J2 = (1 - Ji)/2 (b) J2 = 1 - Ji (c) J2 = 2(1 - Ji) shown 
in Fig. (|3| with parameters ii = 1, L2 = 0, iV = 1024. The 
insets show the corresponding first derivative with respect to 
J I. The diverging discontinuities have exceeded the t/-limits 
of the coordinate axes. 



IV. ENTANGLEMENT 

Given the exact solution of the ID EQCM, we have a 
rare opportunity to analytically probe the entanglement 
in the ground state of a complex system. We here focus 
on one of the most frequently used entanglement mea- 



sure: concurrence [48 ! . The concurrence can quantify 
entanglement for any bipartite system that relates to the 
two-site RMD pij, which is defined as 



C{p^j) = max{0, Ai - A2 - A3 - A4}, (23) 



I^D(P2i-i,2i)(h=1) 
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FIG. 5: (Color online) The derivative of the QD in the Ji — J2 
plane when Li = 1,L2 = 0,A'^ = 1024 for magnetic field 
h — 0.5 and h — 1 respectively. 




FIG. 6: (Color online) The total mutual correlation, the clas- 
sical correlation and quantum discord along the path J2 = Ji 
when Li = 1, L2 = 0, L = 2048 for magnetic field ft = 2. The 
inset shows its first derivative. 



FIG. 7: (Color online) The phase diagram of the extended 
compass model in the transverse magnetic field, Li = 1, L2 = 
0. (I) spin-flop phase for Ji > 0, J2 < 1; (II) antiparallel 
ordered of spin y component for Ji > 0, J2 > 1; (III) sat- 
urated ferromagnetic phase for Ji < 0, J2 > 1; (IV) stripe 
antiferromagnetic phase for Ji < 0, J2 < 1. 



Here p*^ denotes the complex conjugation of pij in the 
standard basis. The eigenvalues of C are 



A 



1,2 



^3,4 



^1 + «) + (a|> + (afa|) 



1 + (<jf 



(25) 

For two spins that are on the same odd bond {2i — 
l,2i}, when Ji > 0, Li = 1, L2 = 0, the ground state 
resides in the subspace, where (crf ) = 0, {(^2i-i'^2i) = 
'o'oi-TO-^i) = (c^i-icrf,), and then it is immediately clear 



'2i-l"2i/ - \"2i 

that C{pij) = — 



2i-l"2i 



On the other hand, when 



Ji < 0, the ground state is in another subspace, where 



[a. 
C{p 
(a 



^3 
2i-l 



) = -{c 



'2il 



0, (a. 



2j/. However, at Ji = 0, 



'2i-l 



2i-l"2i/ 



and then 
0, and then 



C{pij) = 0, which does not support the conjecture in Ref. 
[17] that the first-order QPT at Ji =0 is not signaled by 
the pairwise concurrence. As is displayed in Fig. |8j the 
concurrence of nearest-neighbor pairs of spins captures 
the discontinuity across the first-order transition point 
when there is no external magnetic field, and the diverg- 
ing peaks of the first derivative of concurrence imply the 
continuous QPTs. 



where are the eigenvalues in decreasing order of the 
auxiliary matrix 



(24) 



V. DISCUSSION AND CONCLUSION 

Despite the resemblance in characterizing the QPTs 
in ID EQCM, there still exist much difference between 
the concurrence and the QD. For example, the QD can 
present the correlations between neighbors farther than 
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FIG. 8: (Color online) The concurrence of nearest neighbor 
spins on the odd bonds of with respect to Ji along path 
(a)J2 = 2(1- Ji) (b) Ja = Ji with Li = 1, La = 0, iV = 1024. 
The corresponding first derivatives of concurrence are dis- 
played in insets. 



ture [5U]. The QD even increases with temperature in 
some cases |51j . Furthermore, there is an evidence that 
QD may present a scahng law, which is not the case for 
entanglement General speaking, quantum correla- 

tions are more fundamental than quantum entanglement, 
and may reveal more information about the quantum sys- 
tems. 

In conclusion, we have examined pairwise QD by ex- 
actly solving the ID EQCM in the presence of an exter- 
nal transverse magnetic field. We find that the QD is 
equal to zero for the ID compass model. Remarkably, we 
have successfully extracted information of the location 
and the order of the QPTs in ID EQCM by considera- 
tion of the derivative of the QD with respect to the cou- 
pling parameters. We conclude that a first-order QPT 
is associated with a discontinuity in the QD and a con- 
tinuous second-order transition demonstrates a diverging 
first derivative of QD. The mixed first-order and second- 
order phase transition point features a discontinuity. We 
get the analytic expressions of critical magnetic fields for 
the field-induced QPTs, which are of second order. As a 
result, we then obtain the phase diagram of ID EQCM 
in the transverse magnetic field. For comparisons, we 
show that the pairwise concurrence can characterize the 
phase transitions by exhibiting similar behaviors. Nev- 
ertheless, the QD is believed to be more fundament than 
concurrence in quantifying the quantum correlations, and 
should be of general interest in future studies. 



the next-nearest, while pairwise entanglement may be ab- 
sent for these neighbors {34J. Besides, the QD can charac- 
terize QPTs by exhibiting long-range decay as a function 
of distance in spin systems, which is different from the 
behavior of pairwise entanglement [49\ Moreover, the 
thermal fluctuations extinguish the entanglement, while 
the QD is robust to spotlight QCPs at finite tempera- 
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